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® Process for separating upgoing and downgoing events on vertical seismic profiles. 

@ A process for separating upgoing (45.47) and downgoing (41.43) seismic wave events in a vertical seismic 
profile and particularly for vertical seismic profiles wherein the seismic energy source is offset from the borehole 
in which seismic detectors are located in vertical spaced array. The process operates on two detector signals at 
a time (as Si and Sa) and is based on the concept that waves traveling in opposite directions have spatial 
derivatives of opposite sign. The derivative is approximated by the difference of the two signals which is time 
integrated to recover the phase. The resulting integrated difference signal I is then amplitude scale corrected 
and combined by addition or subtraction with a signal So representing the sum of the two detector signals to 
form a succession of filtered signals which, when recorded in alignment in order of detector depths to form a 
vertical seismic profile preserves either the upgoing or the downgoing seismic events. 
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PROCESS FOR SEPARATING UPGOING AND DOWNGOINQ EVENTS ON VERTICAL SEISMIC PROFILES 



This invention relates to a method of seismic surveying by the vertical seismic profile technique in 
which seismic wave detectors located in a well bore at different depths detect acoustic waves emanating 
from a seismic energy source located at or near the earth surface and generate representative electrical 
signals in response thereto which are recorded in aligned array in order of detector depths, and more 
5 particularly to a method for the processing of seismic data signals for separating the downgoing wave 
events from the upgoing wave events on a vertical seismic profile. 

Seismic sun/eying typically involves the use of a source of seismic energy and an array of seismic 
detectors strategically positioned for its reception. The source of seismic energy can be an apparatus 
capable of delivering a series of impacts or mechanical vibrations to the surface of the earth, the detonation 
10 of a high explosive charge near the earth surface, or other means capable of generating seismic wave 
energy The resulting acoustic waves generated in the earth, including those which are reflected from earth 
strata interfaces, are detected by seismic detectors which transduce the acoustic waves into representative 
electrical signals. From these electrical signals, informational data may be deduced concerning the structure 

of earth's substrata. ^ ^ . ^ u , u 

IS In vertical seismic profiling, the seismic detectors are positioned at different depths in a borehole, such 
as a well bore, and the signals from the detectors in response to reception of energy from a seismic energy 
source are recorded and grouped in alignment in a single display in the order of detector depth. From this 
display, coherences between the signal traces may be noted which may be analyzed and interpreted to 
provide information regarding the geologic substructure. However, because of the multiplicity of compo- 

20 nents comprised in the seismic energy received by a detector and its representative electrical signal, 
analysis of the signals is oftentimes exceedingly difficult. Such components will typically include a 
downgoing component representing the directly arriving wave field propagated from the energy source, 
other downgoing components from seismic energy which have undergone multiple relfections in geologic 
strata above the detector, upgoing components from seismic energy reflected from interfaces of geologic 

25 strata or structures located below the detector, and spurious waves of various kinds. 

Heretofore, a number of field-operating and data processing techniques have been devised to accen- 
tuate the upgoing components in the detector signals and at the same time minimize the interfering effects 

. of the downgoing components since it is the upgoing components and their transit times representative of 
reflections from substrata interfaces which provide most useful information and are of primary concern. 

"0 Most of the current techniques for separating the opposite directi'ons of seismic wave travel in the detector 
signals assume that time shifts which align the first breaks in the signal traces can also be used to align the 
desired upgoing seismic wave events. However, this assumption is strictly true only where the upgoing 
events are derived from horizontal reflectors and the detector signals are from detectors located in a vertical 
borehole with zero source offset, i.e.. the seismic source is located as close to the borehole as possible 

35 rather than being offset therefrom. In particular, the current techniques are of limited effectiveness in 
compensating for "dipping" reflectors and in offset vertical seismic profiling used to identify geologic 
substructures, such as faults or salt donnes located a distance from the well bore. Also, the multichannel 
"dip" filters in current use generally require many more than two detector signals from which to extract 
desired information. 

40 The new filtering process described herein does not require preprocessing to align siesmic events in 
'eismic signal traces and is applicable to a wide range of geometries, including vertical seismic profiling 
applications wherein the seismic source is offset from the well bore. It also makes possible, by means of 
two-trace filtering, to reduce the number of receiver depth levels required in the acquisition of a vertical 
seismic profile, and therefore a reduction in acquisition costs. 

45 Viewed from one aspect the invention provides a method of seismic geophysical surveying wherein a 
seismic energy source positioned at or near the earth surface is adapted to generate seismic wave energy 
v/hich is detectable by seismic detectors in vertical spaced array in a borehole and the electrical signals 
generated by the seismic detectors in response thereto are processed for discriminating the seismic events 
which are upward traveling as they arrive at the detectors from those which are downward traveling as they 

so arrive at the detectors, said method comprising the steps of: 

a) generating seismic waves in the earth at a given location, 

b) detecting the resultant seismic waves at a plurality of vertically spaced levels in a borehole 
penetrating earth formations below the energy source and generating representative electrical signals in 
response thereto by means of detector transducers positioned in the borehole at each said level. 
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c) combining a pair of detector signals from a pair of detectors at adjacent levels in the borehole to 
form a sum signal, 

d) determing the difference between the detector signals to form a difference signal, 

e) integrating the difference signal to form an integral signal, 

5 f) applying an estimated amplitude scale correction to the integral signal to form a resulting signal, 

and 

g) combining the amplitude corrected signal with said sum signal to form a modified filtered signal 
wherein the amplitudes of the seismic wave events which are upwardly traveling as they arrive at the 
detectors are Increased and the seismic wave events which are downwardly traveling as they arrive at the 
10 detectors are attenuated. 

Viewed from another aspect the invention provides a method of seismic geophysical surveying wherein 
a seismic energy source positioned at or near the earth surface is adapted to generate seismic wave 
energy which is detectable by seismic detectors. in vertical spaced array in a borehole and the electrical 
15 signals generated by the seismic detectors in response thereto are processed for discriminating the seismic 
events which are upward traveling as they am've at the detectors from those which are downward traveling 
as they arrive at the detectors, said method comprising the steps of: 

a) generating seismic waves in the earth at a given location, 

b) detecting the resultant seismic waves at a plurality of vertically spaced levels in a borehole 
20 penetrating earth formations below the energy source and generating representative electrical signals in 

response thereto by means of detector transducers positioned in the borehole at each said level. 

c) combining a pair of detector signals from a pair of detectors at adjacent levels in the borehole to 
form a sum signal, 

d) determining the difference between the detector signals to form a difference signal, 
25 e) integrating the difference signal to form an integral signal, 

f) applying an estimated amplitude scale correction to the integral signal to form a resulting signal, 

and 

g) combining the amplitude con-ected signal with said sum signal to form a filtered signal wherein the 
seismic wave events which are downwardly traveling as they arrive at the detectors are preserved. 
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Viewed from a further aspect the invention provides a method of seismic geophysical exploration by the 
vertical seismic profile technique wherein the seismic energy generated by a seismic energy source is 
detected by seismic detectors in vertically spaced array in a borehole and the representative electrical 
signals generated by the detectors in response thereto are recorded in aligned array In the order of detector 
35 depth to form a vertical seismic profile wherein seismic events which are upwardly traveling as they arrive 
at the detectors are separated from downwardly traveling events as they arrive at the detectors, said 
method comprising the steps of: 

a) combining a pair of the recorded detector signals from a pair of detectors at adjacent levels in the 

borehole to form a sum signal, 
40 b) determing the difference between the detector signals to form a difference signal. 

c) integrating the difference signal to form an integral signal, 

d) calculating an estimated amplitude scale correction for application to the integral signal to form a 
resulting signal 

e) combining the amplitude corrected signal with said sum signal to form a filtered signal wherein the 
45 upwardly traveling seismic wave events arriving at the detectors are discriminated from the downwardly 

traveling seismic wave events, 

f) repeating each of the foregoing steps to successive pairs of signals for all the signals of the vertical 

seismic profile, and 

g) recording the successively formed filtered signals in a vertical seismic profile. 
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A preferred embodiment of the invention will now be described, by way of example only, with reference 
to the accompanying drawings, in which: 

Figure 1 is a diagrammatic illustration of an earth cross-section penetrated by a well bore and the 
wave path geometry of seismic waves emanating from a seismic energy source at the earth surface and 
55 received by a seismic detector located at a specified depth in the well bore; 



3 



EP 0 327 758 A2 



Hgure 2 is a- vertical seismic profile display of detector signals representative of seismic waves 
received over time by seismic detectors positions at successively greater depths in a well bore wherein the 
detector signals are arranged in horizontal alignment in correspondence with the depth of the detectors and 
the differing appearances of directly arriving waves at the detectors and wave reflections from interfaces of 
5 aeoloqic strata located below the detectors are illustrated; 

Figure 3 is a reproduction of signal traces 10 and 11 in the vertical seismic profile of Figure 2 where 
these traces are each displayed in horizontal orientation and aligned one below the other; 

Figure 4 is a display, similar to Figure 3. of sum and difference signal traces representing the sum 
and difference of the pair of signal traces 10 and 11 and a third signal trace, an integral trace, representing 
10 an integration of the difference signal trace; 

Figure 5 is a display of a first signal trace resulting from the substraction of the integral trace from 
the sum trace and a second signal trace resulting from an addition of the integral trace and the sum trace; 

Rgure 6 is a vertical seismic profile derived from the vertical seismic profile of Figure 2 after the 
process of the invention has been applied in successive steps, two traces at a time, to all the signal traces 

js in the vertical seismic profile of Figure 2; 

Rgure 7 is a block diagram of a system according to the invention for the filtering of detector signals 

obtained in the acquisition of a vertical seismic profile: 

Rgure 8 is a diagrammatic illustration of the ray-path geometry of an incident wave approaching the 
axis defined by two seismic detectors positioned at two adjacent levels in a well bore; 
20 Figure 9 is a graph in polar form representing the directional amplitude response, (Re. Qo). for the 

filtering process of the invention as a function of incident angle wherein the directly arriving wave arrives at 
an incident angle of O' with respect to the axis defined by the pair of receivers in Figure 7; 

Rgure 10 is a graph in polar form representing the directional amplitude response. R(fl. 6o). for the 
filtering process of the invention as a function of incident angle wherein the directly arriving wave arrives at 
25 an incident angle of 30 * with respect to the axis defined by the pair of receivers in Figure 7; and 

Rgure 11 is a graph in polar form representing the directional amplitude response, R(9, 9o), for the 
filtering process of the invention as a function of incident angle wherein the directly arriving wave arrives at 
an incident angle of 60 * with respect to the axis defined by the pair of receivers in Figure 7. 

30 Referring to the drawings in greater detail. FIG. 1 illustrates schematically a typical arrangement of 
apparatus for conducting seismic operations for the acquisition of a vertical seismic profile. A drilling rig 11 
is shown in position over a well bore 13 which has been drilled in earth formations 15.16.17. The well bore 
13 may be of any specified depth but typically extends to a depth of several thousand feet. A seismic wave 
transducer or detector 19, sometimes called a geophone. is suspended on a conductor cable 21 in the well 

35 bore 13 at a predetermined depth above the interfaces 22.23 of the earth strata 15.16 and 16,17 
respectively The cable 21 which is wound about a reel 24 mounted on the drilling rig 11 or located 
elsewhere on the earth surface near the well bore 13. allows the transmission of electrical signals produced 
by the detector 19 to surface equipment 25 for recording the processing. Such equipment typically includes 
a means (not shown) for providing a precise indication of the depth of the detector 19. 

40 Also shown in FIG. 1 is a seismic wave source 34 located on the earth surface 33 and offset from the 
•veil bore 13 The seismic wave source is depicted in RG. 1 as a weight-dropper 34 which is capable of 
dropping a weight to deliver a series of impacts to the earth, although it is to be understood, that various 
other seismic wave generating means, such as a vibrator or explosive charge, could substitute therefor. 
In operation, seismic waves are generated by an impact delivered by the weight-dropper 34 on the 

45 earth surface 33. These waves are propagated in the earth formations below and detected by the detector 
19 which is fastened against the wall of the well bore 13 by an anchoring element (not shown) in a manner 
oroviding contact between the detector and wail of the well bore which is suitable for the reception of 
seismic wave energy. Such an anchoring element is shown in U.S. Pat. No. 4.397. 004. In response to the 
detection of seismic wave energy, the geophone 19 generates an electrical signal representative of the 

so seismic waves impinging thereon and this signal is transmitted via the conductor cable 21 to a magnetic 
recorder included in the surface equipment 25. While only a single detector 19 is shown in FIG. 1. it is to 
be understood that a plurality of such detectors may be suspended on the cable 21 in vertically spaced 
orientation in contact with the wail of the well bore and their signals transmitted via multiple conductors of 
the cable 21 to the recording equipment which may be any conventional multitrace recorder suitable for 

55 recording the individual geophone output signals in reproducible form. 

It is to be noted frorri observation in FIG. 1 of the ray paths of seismic energy emanating from the 
=ource 31 that some of the seismic wave energy travels downwardly and directly to the detector 19 as 
indicated by the ray path 41 and is normally the first seismic event to be detected by the detector 19. 
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However, there may be other seismic events which are downwardly traveling as they impinge on the 
detector 19. One such, as indicated by the ray path 43. may be a multiply reflected wave representing a 
portion of the seismic energy which strikes the geologic interface 22 and is reflected thereby to reverberate 
between the interface 22 and the earth surface 33. In addition, some of the seismic wave energy which is 
reflected by the interface 22 and illustrated by the ray path 45 is upwardly traveling as it impinges on the 
detector 19 and is the first upwardly traveling seismic event to be detected by detector 19. 

However, the downwardly traveling seismic wave energy which striked the interface 22 is not totally 
reflected therefrom and some of the seismic wave energy continues downwardly to be reflected by other 
geologic interiaces below the interface 22, such as the interface 23. A second reflection of seismic wave 
energy as indicated by the ray path 47 is upwardly traveling as it impinges on the detector 19 and is the 
second upwardly traveling event to be detected. There may, of course, be several more of such upwardly 
traveling seismic waves representing reflections from lower geologic strata which are not shown in FIG. 1 
for clarity of explanation. ^ ^ ^ 

Hgure 2 is an illustration is a synthetic vertical seismic profile display of receiver depth and travel time 
showing how the seismic events associated with the ray paths in FIG. 1 would appear in a vertical seismic 
profile. Each of the vertically oriented traces in the profile corresponds to the recorded electrical signal 
transmitted by a single detector in the well bore and the total 51 traces are recorded signals from a detector 
at 51 successively lower depths in the well. However, it is also to be understood, that the vertical seismic 
provUe could be acquired with only a single generation of seismic energy from the energy source using 51 
detectors in the well bore or by a series of impacts of the weight dropper wherein one or more detectors 
are successively positioned at greater depths in the well bore to receive the energy. 

In the vertical seismic profile of Rgure 2. the downgoing events on the provile are the directly arnvmg 
wave from the seismic energy source and a multiple reflected wave revertserating in the first layer of earth 
below the energy source whereas the upgoing events are the reflections of seismic energy from two 
interfaces of earth strata. These events appear on a trace as a sudden increase or a peal< in signal 
amplitude and in a vertical seismic profile these amplitude changes in the traces of the profile are disposed 
as a dipping feature of the profile. It will be noted that the amplitude changes representing the downgoing 
arriving waves in the several traces appear in alignments which are dipped oppositely to the alignment of 
the signal amplitude changes representing the arrival of upgoing reflections which phenomenon Is typical of 
vertical seismic profiles wherein opposite directions of seismic wave travel within the earth lead to opposite 
dips on the display. 

Sin ce e x p lorat i on by v o rtical s oi smic profiling techninues i. s concerned with detection of the reflecting 

horizons of the different earth strata, the detection and analysis of the upgoing waves is the primary 
concern. It is therefore necessary, in order to improve the signal-lo-noise ratio and to better distinguish the 
upgoing waves, to filter the collected signals. It is therefore desirable that the filtering prdcess operates to 
suppress the downgoing events in the vertical siesmic profile. 

The filtering process of this invention will be hereinafter described by consideration of traces 10 and 11 
in the vertical seismic profile of FIG. 2. These traces have been selected so that the undesired multiple and 
the second reflection interfere with one another in that they add constructively on the first trace and cancel 
on the second. Despite such interference, it will nevertheless be shown that the process of this invention 
can correctly separate the two events. 

In FIG 3 the traces 10 and 11 of FIG. 2 are reproduced exactly but are plotted in horizontal orientation 
for purposes'of explanation. From observation of FIG. 3 it will be noted that the wavelet shape for the first 
two events on each trace is a trough followed by a peak. The second two events, corresponding to the 
second reflection and the multiple in the vertical seismic profile of FIG. 2. interfere with one another to 
construct a different wave shape (trough-peak-trough) in the first trace (trace 10) and to perfectly cancel one 
another on the second trace (trace 1 1 ). 

The first step in the process is to form two new signals from the sum and difference of the original 
traces 10 and 11. These new signals are represented in FIG. 4 by the signal traces So and (S1-S2). From 
so observation of these traces, the summed or averaged trace So. with respect to wave shape, looks much like 
the original traces 10 and 11. However, the difference trace (8,-82) has a different wavelet shape for the 
events appearing thereon in that each wavelet is approximately 90* phase shifted, and. more important, 
shows a change in polarity between the upgoing and downgoing events. It is therefore to be noted that the 
first two events on the summed trace, one desired and the other undesired, look identical whereas on the 
55 difference trace, with respect to these same two events, one is upside down relative to the otiier. 

The process then requires a wave shape correction which is achieved by the step of integrating the 
difference signal to obtain an integral thereof, represented in FIG. 4 as the integral trace I. An amplitude 
scale correction of the integral signal is then performed followed by the step of substracting the integrated 
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difference signal from the summed signal. The subtracting step in effect doubles the amplitudes of the 
upgoing events and cancels the downgoing events, as is readily apparent in the signal trace (So -I) shown in 
FIG 5. On occasion, the arrival of the downgoing events may be of interest and in such case, the 
downgoing events may be preserved by the step of adding the integrated and summed signals, So and I. to 
5 obtain the resulting signal trace (So + 1) shown in FIG. 5. For either of the traces of FIG. 5, the arrival times 
of the separated events are appropriate for a depth halfway between the receiver pair. Also, the amount of 
residual unwanted events on the filtered signal traces is a function of the receiver spacing and seismic 
wavslet 

The filtering process of this invention is ideally suited to processing vertical seismic profiles, but has 
10 potential application wherever a dip filter is desired. It can be used whenever desired and undesired events 
are dipping in opposite directions, or can be shifted so as to appear dipping in opposite directions. For 
example, a surface wave may be shifted relative to reflection energy in order to suppress the surface wave. 

As previously stated, there will normally be a requirement for applying an amplitude scale con-ection to 
the integral signal trace before performing the step of subtracting the integral trace I from the sum trace So. 
75 This con-ection. of course, can be specified exactly if the time separation, as might be available from an 
acoustic well log or from the vertical seismic profile first breaks, is known. Nevertheless, the scaling can 
effectively be achieved by power balancing the traces before subtraction. 

A procedure whereby a time delay scaling factor can be estimated with appropriate accuracy is herein 
described. 

20 In the acquisition of a vertical seismic profile, let Si and Sa be the signals detected at two adjacent 
levels level 1 and level 2. in a well bore. Also, let the signal which would have been measured at the 
midpoint of the adjacent levels be denoted by So. If T be the time delay of the wave traveling between the 
midpoint and the levels for acquisition of Si and Sz. then, assuming a constant velocity layer. 
So = u(t) + d(t) (1) 

25 Si = u(t - T) + d(t + T) (2) 

$2 = u(t + T) + d(t - T) (3) . . , , 

•Miere u(t) and d(t) are the upgoing and downgoing wave fields at the midpoint level. 

Si and Sz ate then expanded in a Taylor's series wherein, for closely spaced repeiyers. terms of first 
order in T are retained to obtain 

30 S- = u - Tu + d + Td (4) 
Sz = u + Tu + d - Td (5) 
wherein the dot denotes differentiation with respect to t. 

The average of Si and Sa is therefore 
S = (Si + S2)/2 = u + d (6) 

35 The difference of Si and S2 is 
S2 - Si = 2Tu - 2Td (7) 

The integral of the difference is 
I - 2T (u-d) (8) 

In equation (8), the negative coefficient of d is due to the opposite directions of travel of the two wave 
40 fields. The factor 2 comes from the definition of T to be the time delay between the midpoint and level 1 or 
2, that is. the time delay between level 1 and level 2 is 2T. 

Dividing I by 2To, where T, is some arbitrary constant, and adding the result to S. the following is 

obtained: 

Xiizdi + (u+d) 



or ^^S-O u . (1 -^^) d (9) 
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It is therefore apparent that u and d can be determined from the combination of S and I expressed in 
equation (9). The downgoing wave field may be eliminated by setting To = T to obtain u. Alternatively, the 
upgoing wave field may be eliminated by setting To = -T to obtain d. 

T can be determined from the difference in the first arrival times of a wave field at the two detectors, 
since the difference represents the travel time between the two detectors. This travel time divided by 2 
therefore equals T. However, a more flexible and automatic procedure for determining T is outlined below. 
From equations (6), (7) and (8), the power in S and I is given by 

= Pu + Pd + 2 ^.o (u.d) (10) 
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p, - 4T2 [P„ + Pd - 2 *o (u,d)j (11) o- u 

where *o (u-d) is the zero lag cross-correlation of the upgoing and downgo(ng wave fields. Since the 
upgoing and downgoing wave fields at a given time are uncorrelated. «o(u.d) = 0, Then, by dividing 
equation (11) by equation (10). the following is obtained: 

5 4 T2 = P|/Ps (12) 

Accordingly T which is readily obtainable from equation (12). provides an estimated value of the time delay 
between the two levels, which estimate is not based on a single value, but represents an average value of 
the shifts for the whole trace. This estimated value is therefore a better estimate than can be determined 
from the first arrivals alone. , ^. „ u 

w It is to be noted, however, that this derived estimated value of T is valid so long as T is small enough 
that terms higher than first order in the Taylor's series can be neglected. A small enough T means that the 
receivers are sufficiently close together that T is small compared to the smallest period (associated with the 
highest frequency) present in the data. If T becomes a large fraction of the smallest period, then a condition 
similar to "cycle skipping" occurs. That is. when comparing two cosine waves which have a relative time 

IS shift it is impossible to distinguish between waves shifted up greater than half a cycle and down less than 
half a cycle. Since the process of this invention attempts to detect these relative shifts, the process will fail 

for those wavelengths. . ^ . ,^ 

From the foregoing, it is to be concluded that the filtering process herein descnbed yields a resulting 
signal approximating the signal which would have been recorded midway between the two receivers, in 

20 effect the filter process mixes two traces whereas most multichannel filters mix several traces. This, of 
course results in the loss of one trace per group of input traces and also results in a new depth 
specification. However, since there is normally less level smearing with a two-trace filter, it better preserves 
faults and steeply dipping events in offset vertical seismic profiling. 

Figure 6 shows the results of applying the process of the invention two-by-two to all the traces in the 

25 synthetic vertical seismic profile of FIG. 2. It is shown therein that a single pass or application of the 
process greatly enhances the upgoing events relative to the downgoing events. In effect, the upgoing 
events are doubled in amplitude and the downgoing events are so attenuated that only a small residual 
could be noted. A second pass would further improve the separation but normally is not needed unless 
significant differences in dip occur. For example, if mode converted events are present, one pass would 

30 reject the downgoing P-wave (compressional wave) energy and a second pass might be needed to reject 
the downgoing S-wave (shear wave) energy. ...^ , ^ . . 

A system for implementation of the filtering method described above is shown in FIG. 7. The detector 
signals obtained from vertically spaced detectors in a borehole are recorded and stored in a memory 60. 
The memory 60 is controlled to furnish for processing a pair of detector signals. S. and Sz. corresponding 

35 to signals from detectors at adjacent levels. The pair of signals are also delivered to an adder 61 and a 
subtracter 62 for obtaining sum and difference signals. So and (Si-Sj). respectively. The difference signal is • 
then applied to an integrating circuit 63 which produces an integral signal I. The integral signal I is then 
delivered to a calculator 65 which applies an amplitude scale correction to the integral signal I as 
determined from the equation (12) above. 

40 The amplitude-corrected integral signal and the sum signal are delivered to a subtracter circuit 67 
wherein the amplitude-corrected integral signal is substracted from the sum signal So to provide an output 
signal So - l.c wherein the amplitudes of upgoing seismic events are significantly enhanced and the 
downgoing events are significantly attenuated. The output signal from the subtracter 67 is then recorded in 
memory 68 from which it may be applied to a recorder display device 69 for display in a vertical seismic 

'''°Tftematively. or in conjunction therewith, the sum signal So and the amplitude-corrected integral signal 

may also be delivered to an adder 71 wherein they are added to presen/e the downgoing events, and the 

resulting summation signal stored in a memory 72. 

It is also to be understood that the system is adapted to process all the signals in a vertical seismic 
50 profile such that the process is repeatedly applied two-by-two to all the signal traces of the profile to obtain 

an improved vertical seismic profile in accordance with the invention. The system can be used with either 

analog or digital means of data processing. 

An important consideration is the response of the filter to more than one dip. Accordingly, it will herein 

be shown that the time scale factor which is chosen to eliminate a particular downgoing event also 
ss attenuates similar downgoing events and enhances all upgoing events. In this respect, it is to be noted that 

the scaling factor T can be related to the angle of approach of the incident wave. Referring to FIG. 8, which 

diagrammatically illustrates the approach to a receiver of an incident wave arriving at an angle to the axis 

line defined by a pair of aligned receivers, it will be seen that 
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2T = AZC0S(9)/V (13) . , , 

whereAZ is the distance between two levels. V is the velocity of the wave in the region of the two levels, 
and 9 is the angle between the direction of the wave travel and the axis defined by the receiver pair at two 
adjacent levels. A wave traveling straight down the axis corresponds to 9 = 0 . 
5 ■ Assuming a scaling factor 2 To is chosen to reject a downgoing wave of such that 2 1, =AZ cos 
(9o)/V. then if the upgoing wave has 2 T = AZ cos (9)/V. the response curve for the filter as obtained from 
equation (9) is given by 

R (9. «o) = 1 - T/To = 1-cos(e)/cos(ao) (14) . . n ^ • c- ..n 

The graph of equation (14) is shown in Rgure 9 for an incident wave where do = 0 and in Figures 10 
10 and 11 where So = 30* and 60* . respectively. In these curves, the distance from the origin corresponds to 
H(e 9o) and the beam pattern polar plots indicate amplitude by circles of successively greater radius as a 
function of incident angle indicated by radial lines. In each of the plots, the vertical axis corresponds to well 
bore geophone alignment and the arrows indicate wave direction used for time scaling, in each of Figures 9. 
10 and 11. the curve for R(s. 9,) traces the filter, response for seismic wave events arriving from ail other 

IS directions. . . ^ .i. ■ • * 

It will be noted from Figure 9. wherein ^0= 0 . that a wave approaching the axis of the receiver pair at 
an angle of 0* as represented by the upper arrow is cancelled exactly, with the response curve depicting 
an amplitude of zero. Also, it is to be noted that a wave, such as a reflected wave, which approaches the 
receiver pair at an angle of 180* as represented by the lower arrow is enhanced by a factor of 2. Similarly. 

20 a wave event which approaches the receiver pair at an angle of 90*. i.e. broadside to the axis of the 
receiver pair, arrives at both receivers at the same time and is passed exactly, the amplitude being 
represented equal to 1 . 

Referring to Figures 10 and 11. the upper arrows might typically represent the direct arrivals of seismic 
waves from an offset seismic wave source and the lower arrows, a reflected wave from a horizontal 
25 substratum It will be observed that the direction of exact cancellation is such that it is always at and the 
magnitude of the response curve for a wave arriving at 180* . R(180* - 9o.0o). is always 2. From observation 
of each of the curves of FIGS. 9, 10 and 11, it will be noted that waves which approach the receiver array 
from below are always enhanced (R>1). while waves which approach the array in a wide band around 0 
are always reduced (R<1). 

30 It is to be noted, therefore, that the two-trace filter of this invention has a good rejection response over a 
wide range of dip angles. It is particularly effective when the number of traces is small and for dip moveouts 
that are not linear, as in offset vertical seismic profiling. However, if several different dips are to be rejected, 
a second pass of the filter will prove advantageous. 

Applications have shown the two-trace directional filtering process of the invention to be tolerant of poor 

35 «ignal-to-noise ratios and will not produce artificial alignments like some harshly applied filters. However, 
since the process will enhance low frequencies relative to highs, it is generally advisable that a band-pass 
filter and/or spectral balancing operator be subsequently applied. Also, since in the processing of vertical 
seismic profiles, the invention does not require first break time picks to accomplish separation of the 
upgoing and downgoing events, the process can be applied at an earlier stage than other dip filters. 

40 It is not a requirement of the invention that receiver spacing must be uniform, although the spacing 
must be small enough not to alias frequencies of interest. If the spacing selected is too large compared to a 
wavelength of seismic energy, then wavelet distortion will result which is true, in general, for all dip filters. 
While the filter process of the present invention is useful for any wavelengths which are unaliased. the filter 
process is preferably used with about half the nomnal range of spacing, i.e., best results are achieved with 

45 spacings less than a quarter of the shortest wavelength. Successful applications of the process have been 
achieved with vertical seismic profiles wherein the level spacings are between 10 and 125 feet. Accordingly, 
the technique can be used effectively on the majority of existing vertical seismic profile data sets. 

The technique of the invention also allows for large depth separations between closely spaced pairs of 
receiver locations and consequently a significant reduction in the total number of detector levels as 

=0 rompared to most existing techniques which require a small depth spacing between all detector levels and 
oftentimes produce more levels than are needed. The larger number of detecting levels are normally 
associated with greater costs in service charges and rig time. Furthermore, since it is possible to practice 
the invention using a ten foot detector level spacing, it also becomes possible to place the receiver pair m 
the same tool, thereby allowing even more acquisition flexibility and further possibilities in reducing costs. 

55 From the foregoing description of the invention, it is to be appreciated that the two-trace filter process 
of this invention has far more general applications than existing techniques. There are many potential 
advantages of the two-trace filter which lie in the less restrictive assumptions imposed on acquisition 
geometry and structural complexity of the subsurface geology, such that potential applications are to be 
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found wherever a dip filter is desired. 

It is also to be understood that alternative techniques in performing the various method steps of the 
invention may be used as are apparent to those skilled in the art. without departing from the scope of the 
invention as defined in the appended claims. 



Claims 

1 A method of seismic geophysical surveying wherein a seismic energy source positioned at or near 
10 the earth surface is adapted to generate seismic wave energy which is detectable by seismic detectors in 
vertical spaced array in a borehole and the electrical signals generated by the seismic detectors in 
response thereto are processed for discriminating the seismic events which are upward traveling as they 
arrive at the detectors from those which are downward traveling as they arrive at the detectors, said method 
comprising the steps of: 
1$ a) generating seismic waves in the earth at a given location, 

b) detecting the resultant seismic waves at a plurality of vertically spaced levels in a borehole 
penetrating earth formations below the energy source and generating representative electrical signals in 
response thereto by means of detector transducers positioned in the borehole at each said level. 

c) combining a pair of detector signals from a pair of detectors at adjacent levels in the borehole to 

20 form a sum signal. 

d) determining the difference between the detector signals to form a difference signal. 

e) integrating the difference signal to form an Integral signal, 

f) applying an estimated amplitude scale correction to the integral signal to form a resulting signal. 

25 ^'^ g) combining the amplitude corrected signal with said sum signal to form a modified filtered signal 
wherein the amplitudes of the seismic wave events which are upwardly traveling as they arnve at the 
detectors are Increased and the seismic wave events which are downwardly traveling as they arrive at the 
■ detectors are attenuated. 

30 2 A method as recited in claim 1 wherein the seismic energy source is offset from the borehole. 

3 A method as in claims 1 or 2 wherein the electrical signals from the plurality of detectors are first 
recorded In aligned array in order of detector depth to form a vertical seismic profile prior to the step of 
combining a pair of detector signals. .... 

4 A method as recited in claims 1. 2 or 3 wherein the seismic wave energy generated by the seismic 
35 energy source is detected by vertically spaced pairs of seismic detectors located at successively lower 

depths in the borehole and the steps of the method are successively applied to the signals of each said 

^^"^5 A method as recited in claim 4 wherein the filtered signals from the processed pairs of detector 
signals are recorded in the order of their successively greater depths to form a filtered vertical seismic 

40 ProWe.^ ^^^^^ geophysical sun/eying wherein a seismic energy source positioned at or near 

the earth surface is adapted to generate seismic wave energy which is detectable by seismic detectors in 
vertical spaced array in a borehole and the electrical signals generated by the seismic detectors in 
response thereto are processed for discriminating the seismic events which are upward traveling as they 

45 arrive at the detectors from those which are downward traveling as they arrive at the detectors, said method 
comprising the steps of: 

a) generating seismic waves in the earth at a given location, 

b) detecting the resultant seismic waves at a plurality of vertically spaced levels in a borehole 
penetrating earth formations below the energy source and generating representative electrical signals in 

50 response thereto by means of detector transducers positioned in the borehole at each said level, 

c) combining a pair of detector signals from a pair of detectors at adjacent levels in the borehole to 
form a sum signal, 

d) determining the difference between the detector signals to form a difference signal. 

e) integrating the difference signal to form an integral signal, 

f) applying an estimated amplitude scale correction to the integral signal to form a resulting signal, 

and 



55 



9 



EP 0 327 7S8 A2 



10 



IS 



20 



25 



30 



35 



40 



45 



g) combining tlie amplitude corrected signal with said sum signal to form a filtered signal wherein the 
seismic wave events which are downwardly traveling as they arrive at the detectors are preserved. 

7 A method as recited in claim 6 wherein the seismic energy source is offset from the borehole. 

8 A method as recited in claims 6 or 7 wherein the seismic wave energy generated by the seismic 
energy source is detected by vertically spaced pairs of seismic detectors located at successively lower 
depths in the borehole and the steps of the method are successively applied to the signals of each said 

''^"^9 A method as recited in claim 8 wherein the filtered signals from said successive pairs of detector 
signals are recorded in the order of their successively greater depths to form a filtered vertical seismic 

"^"^^'lO In a method of seismic geophysical exploration by the vertical seismic profile technique wherein the 
seismic energy generated by a seismic energy source is detected by seismic detectors in vertically spaced 
array in a borehole and the representative electrical signals generated by the detectors in response thereto 
are recorded in aligned array in the order of detector depth to form a vertical seismic profile wherein 
seismic events which are upwardly traveling as they arrive at the detectors are separated from downwardly 
traveling events as they arrive at the detectors, said method comprising the steps of: 

a) combining a pair of the recorded detector signals from a pair of detectors at adjacent levels in the 

borehole to form a sum signal. 

b) determining the difference between the detector signals to form a difference signal. 

c) integrating the difference signal to form an integral signal. 

d) calculating an estimated amplitude scale correction for application to the integral signal to form a 

resulting signal, . , ^ . ^. 

e) combining the amplitude corrected signal with said sum signal to form a filtered signal wherein the 
upwardly traveling seismic wave events arriving at the detectors are discriminated from the downwardly 

traveling seismic wave events. ..... , ^ ^- . 

f) repeating each of the foregoing steps to successive pairs of signals for. all the signals of the vertical 

seismic profile, and 

g) recording the successively formed filtered signals in a vertical seismic profile. 
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